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Subsonic propellers in a strong wind as anomalous X-ray pulsars 
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Abstract. The appearance of subsonic propellers situated in a strong wind is discussed. We show that it is similar to the 
appearance of anomalous X-ray pulsars (AXPs) provided the mass and the magnetic moment of neutron stars are M = 1.4Mq, 
and/i = 2 x 10 30 Gem 3 , respectively, and the strength of the wind is M c = 3x 10 17 g s -1 . Under these conditions, the spin periods 
of subsonic propellers are limited within the range of 5-15 s, and the expected spin-down rates are close to P = 7 x 10~" s s -1 . 
The mass accretion rate onto the stellar surface is limited to the rate of plasma penetration into its magnetosphere at the 
boundary. As this process is governed by the reconnection of the field lines, the accretion rate onto the stellar surface constitutes 
1-2% of M c . In this case the X-ray luminosity of the objects under consideration can be evaluated as Lx ~ 4 x 10 35 erg s . The 
model predicts the existence of at least two spatially separated sources of the X-ray emission: hot spots at the stellar surface, 
and the hot atmosphere surrounding the magnetosphere of the star. The ages of the subsonic propellers under the conditions of 
interest are limited to S> 10 s yr. 
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1. Introduction 

Anomalous X-ray pulsars (AXPs) are currently recognized as 
a separate class of single magnetized neutron stars, with com- 
mon characteristics different from those of the normal pulsars. 
The theoretical interpretation of these objects is still a point of 
discussion. A successful theoretical model of AXPs should be 
able to explain the following basic properties of these objects 
(see Mereghetti et al. 2002 for a comprehensive review): 

(1) the clustering of spin periods in the range 5-12 s; 

(2) the spin-down rates of 10~ 10 - 10~ 12 s s _1 ; 

(3) the soft X-ray spectra, which are clearly different from those 
of the accretion-powered pulsars; 

(4) the X-ray luminosity at a level of ~ 10 35 - 10 35 erg s _1 ; 

(5) relatively young ages of these objects, derived from the as- 
sociation of at least two of AXPs with supernova shell remnants 
(SNRs). 

A key question in the modeling of AXPs is the state of 
the neutron star. Actually, there are only 4 possible answers: 
ejector, supersonic propeller, subsonic propeller, and accre- 
tor (see e.g. Ikhsanov 2001a). Three of these possibilities have 
already been discussed in the literature. The state of ejector (= 
spin-powered pulsar) was analyzed within the magne tar model. 
The observed spin-down rates, the soft spectra and the young 
ages of AXPs can be interpreted within this model provided the 
strength of the dipole magnetic field of the star is B ^ 10 15 G 
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(see e.g. Thompson & Duncan 1996 and references therein). 
At the same time, the clustering of the periods can hardly be 
explained in terms of magnetars unless some deviations from 
the canonical scheme of the evolution of neutron stars are in- 
voked in the model (for discussion see Psaltis & Miller 2002 
and Colpi et al. EBDOl . 

The excess of the observed luminosity over the spin-down 
power allows us to reject the possibility that AXPs are super- 
sonic propellers, i.e. spinning-down neutron stars for which 
the condition R COI 55 R m < min{R a ,R\ c } is satisfied. Here 
R COI = (GM n jLo 2 y /3 , R a = 2GM ns /y^ 1 and R m are the corota- 
tion, accretion, and magnetospheric radii of the neutron star, 
respectively. R\ c - c/w is the radius of the light cylinder. 
ai = 27r/P s and M ns denote the angular velocity and the mass 
of the star, and V re i is the relative velocity between the star and 
the surrounding medium. G and c are the gravitational constant 
and the speed of light. On this basis, the corresponding state is 
presently discussed (e.g. Marsden et al. 2001 1 only in the con- 
text of a previous spin-down evolution of AXPs to their present 
state. 

The modeling of AXPs as neutron stars accreting mate- 
rial from a disk has been reported by many authors (see e.g. 
Chatterjee, Hernquist and Narayan 2000 and Mereghetti et al. 
2002 for a detailed description). Within this approach, the ob- 
served luminosity is interpreted in terms of the accretion power. 
Therefore, the 'power deficit' problem mentioned above turns 
out to be naturally avoided. However, to meet other criteria, 
such as the relatively high spin-down rates, the soft X-ray spec- 
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trum and the clustering of periods in a relatively narrow range, 
a rather specific scheme of neutron star evolution has to be in- 
voked (Eksi & Alpar l2"UU3l . 

In the light of this situation, a consideration of the last pos- 
sibility, namely, that the state of neutron stars in AXPs is sub- 
sonic propeller appears quite reasonable. The corresponding 
analysis is the subject of this paper. As shown below, the mod- 
eling of AXPs within the subsonic propeller approach appears 
effective in several important aspects. In particular, it provides 
us with an explanation of the observed clustering of periods 
(Sect. 2) and spin-down rates (Sect. 3). It also predicts the ex- 
istence of two spatially separated sources of X-ray emission 
(Sect. 4) with the total luminosity of order of the observed lu- 
minosity of AXPs. The ages of supersonic propellers are eval- 
uated in Sect. 5. The proposed scenario is discussed in Sect. 6 
and the basic conclusions are summarized in Sect. 7. 

2. Period clustering 

As recently shown by Ikhsanov ( 2001b I, the spin period of a 
neutron star in the state of subsonic propeller satisfies the fol- 
lowing condition: P c d 5= P s SS Pbr, where 

P o6 c*3.2(j.% 1 m- 5 ' 7 M^ n s, (1) 

is the period at which the star switches its state from supersonic 
to subsonic propeller, and 

P bl - lep^M-^m-^s, (2) 

is the period at which the state of the star switches to accre- 
tor. Here m and p^ are the mass and the dipole magnetic mo- 
ment of the star expressed in units of Mq, and 10 30 Gcm 3 , re- 
spectively. Mn is the strength of the wind, M c , expressed in 
units of 10 17 g s _1 . The meaning of this parameter is the mass 
of the surrounding material interacting with the neutron star, 
moving through the medium with the velocity V re i, in a time 
unit: M c = nR^poaV^x. Here p m is the density of material sur- 
rounding the star at the distance R a . 

Combining Eqs. Q and © we find that the clustering of 
periods of AXPs can be explained in terms of the subsonic pro- 
peller approach provided these objects are neutron stars of mass 
M ~ 1.4 Mo and magnetic moment p ~ 2 x 10 3() Gcm 3 , situ- 
ated in a wind of strength M c ~ 3 x 10 17 g s _1 . 

3. Spin-down rate 

As shown by Davies & Pringle ( 198 1 1, the spin-down power of 
stars in the state of subsonic propeller can be evaluated as 

L ssp = 8 x 10 33 /4 in 1 P- 3 egr s" 1 , (3) 

where P5 = P s /5 s. This means that subsonic propellers are 
expected to spin down at a rate 

P 3 L 

P = -^j- ~ 2.5 x 10-' 1 (4 m~ x I A l s s -1 , (4) 

where 745 is the moment of inertia of the neutron star expressed 
in units of 10 45 g cm 2 . 




Fig. 1. Schematic description of the subsonic propeller. A mag- 
netized neutron star rotates with angular velocity co. The mag- 
netosphere of the star is surrounded by a hot turbulent enve- 
lope. As the star moves through the medium with relative ve- 
locity V re i, the gas overflows the outer edge of the envelope 
compressing the envelope plasma. The overflow rate is equal 
to the wind strength, M c . The rate of plasma motion through 
the envelope is equal to the rate of plasma penetration into the 
magnetosphere of the star at the inner boundary of the enve- 
lope, M lec . The material penetrating into the magnetic field of 
the star flows along the field lines and reaches the stellar surface 
at the magnetic poles. The hatched regions represent shocks in 
front of and behind the subsonic propeller moving through the 
dense environment. 

Putting the above derived parameters (M ~ 1 .4 Mq and p. ~ 
2 X 10 30 G cm 3 ) into Eq. 10}, one finds that the spin-down rate 
of AXPs within the subsonic propeller approach is expected to 
be of order 7 X 10 -11 s s" 1 . This is in a good agreement with the 
typically observed values. 

4. Accretion power 

According to the picture presented by Davies & Pringle ( 198 1 1, 
the magnetosphere of a neutron star in the state of subsonic 
propeller is surrounded by an adiabatic (p oc R~ 5 ^ 2 ) spherically 
symmetrical plasma envelope. The energy input to the enve- 
lope due to the propeller action by the star dominates the en- 
ergy losses from the envelope as long as P s < P^. Therefore, 
the temperature of the envelope plasma remains of order of 
the free-fall temperature, T p i - T$ = GM ns m p /kR and, corre- 
spondingly, the sound speed is of order of the free-fall velocity, 
V s - Vg = y/2GM 1K /R (here m p , and k are the proton mass and 
the Boltzmann constant). Under these conditions the height of 
the homogeneous atmosphere throughout the envelope is com- 
parable to the radius, R, and the envelope is extended from the 
magnetospheric boundary up to the accretion radius of the star. 

Within this picture (see Fig.0 the envelope is quasi-static: 
the characteristic time of the drift of material through the enve- 
lope in the radial direction is much larger than the characteristic 
time of the turbulent motions. The rate of the drift is signifi- 
cantly smaller than M c , i.e. the rate by which the circumstellar 
material overflows the outer edge of the envelope compressing 
the envelope plasma (for discussion see Ikhsanov 2003, and 
references therein). 

Under the condition P s > P c d, the centrifugal barrier is not 
effective: the plasma, penetrating from the envelope into the 
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magnetic field of the star, is able to flow along the magnetic 
field lines and to accrete onto the stellar surface. Therefore, 
an accretion process onto the stellar surface can be expected. 
However, the rate of this process is limited to the rate of plasma 
penetration into the magnetic field of the star. 

As shown by Arons & Lea (JT976a I and Eisner & Lamb 
J1976i . the interchange (e.g. Rayleigh-Taylor) instabilities of 
the magnetospheric boundary under the condition T p i ~ Tff are 
suppressed (we would like to note that the observed luminosi- 
ties of AXPs are a factor of 5-10 lower than the critical lumi- 
nosity for which Compton cooling is effective). In this situa- 
tion, the process of plasma penetration into the magnetosphere 
is governed by the Bohm diffusion and/or the reconnection of 
the field lines at the magnetospheric boundary. An evaluation 
of the corresponding rates of penetration has recently been re- 
ported by Ikhsanov (2001a, see Eqs. 22 and 24). Using these 
results, one finds that the rate of plasma penetration into the 
magnetosphere due to the reconnection of the field lines is high 
enough for the X-ray luminosity of subsonic propellers to be of 
the order of that observed from AXPs (see Eq. (24) in Ikhsanov 
I2()()lat 



^acc — "**re 



GM n , 

fins 



4 x 10 35 ergs 



(5) 



1.4M© 



10 6 . 



0.1 



0AR„ 



M, 



3xl0 17 gs- 



Here M rec is the rate of plasma penetration from the envelope 
into the magnetosphere due to the reconnection process, and 
/? ns is the radius of the neutron star. a r is the efficiency of the 
reconnection process, which is normalized according to the re- 
sults of modeling of the reconnection processes in solar flares 
and in the Earth's magnetopause (see Priest & Forbes 2000|and 
references therein). A m is the average scale of plasma inhomo- 
geneities at the base of the envelope, which is normalized fol- 
lowing Arons & Lea i 1976 b). and Wang & Robertson I1985> . 

This indicates that the radiation of subsonic propellers sit- 
uated in a strong wind is powered mainly by the reconnection- 
driven accretion of material onto the stellar surface. At the same 
time, from 5 to 10 percent of the radiation emitted by these ob- 
jects are expected to be contributed by the hot envelope sur- 
rounding the magnetosphere (see Eq.[3}- Therefore, within the 
considered approach the existence of at least two spatially sep- 
arated sources of the X-ray emission can be expected. 

A detailed analysis of the radiation spectrum of subsonic 
propellers is beyond the scope of this paper and will be pre- 
sented in a forthcoming paper. 



5. Age 

The age of subsonic propellers can be expressed as the super- 
position of pulsar-like and supersonic propeller spin-down time 
scales. 



5.1. Pulsar-like spin-down time scale 

The period at which the neutron star switches its state from 
ejector to supersonic propeller is (see Eq. 6 in Ikhsanov 2001ai 



P mi = 0.2 /i^ M^ 4 V- 1 ' 4 s. 



The spin-down rate of the star in the ejector state is 
87rV 



P = 



3c 3 IP 



(6) 



(7) 



Combining these equations, one finds the spin-down time scale 
(T a = P/2P) as 



r a * 1.2 X 10 5 ^.3 7 45 M n '/ 5 2 V- l/2 yr, 
where Vj = y re i/10 7 cms _1 . 



(8) 



5.2. Supersonic propeller spin-down time scale 

The spin-down rate of a star in the state of supersonic propeller 
is 

P 3 L 



4n 2 I 



(9) 



where L sups is the rate of the rotational energy loss by the star in 
the state of supersonic propeller. In the general case, this value 
can be limited to $ MGM/R m . Therefore, combining Eqs. 
and (j9Ji, one can express the lower limit to the spin-down time 
scale of stars in the supersonic propeller state as 



t c - 8 x IOV303 '45 M-^^yr. 



-11/14 t/1/2 , 



(10) 



Although our estimate represents the absolute minimum for t c , 
the real value of this parameter is unlikely to be significantly 
larger. In particular, as shown by Lipunov & Popov ( 1995 1, the 
duration of the supersonic propeller state under all reasonable 
conditions turns out to be smaller than the spin-down time scale 
of the star in the state of ejector. On this basis, one can limit the 
age of the star as it has evolved to the state of subsonic propeller 
as 



(Ta + T c ) < T ^ 2T a 



5.3. Subsonic propeller spin-down time scale 



(11) 



Finally, the spin-down time scale of the star in the subsonic 
propeller state can be expressed following Ikhsanov (2001b 
Eqs. 10), as 



fd-2x 10 3 // 3 o 3 7 45 P5 m.4 yr, 



(12) 



where mi 4 - \ AMq. 

Ta is small compared to the value of both, r a and t c . Hence, 
the age of the neutron stars with the parameters derived in 
Sect. 2 is determined mainly by the pulsar-like spin-down time 
scale and can be limited to < 2r a . In order to fit the currently 
adopted evaluations of the ages of AXPs, one has to suggest 
that in a previous epoch the conditions in the circumstellar en- 
vironment were slightly different. In particular, one could as- 
sume that the strength of the wind in the previous epoch was 
about 10 18 g s _1 , which corresponds to the accretion process at 
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the Eddington limit. On the other hand, one can also envisage 
a situation in which the value of the relative velocity V K \ is 
slightly larger than 100km s _1 (this normalization of the ve- 
locity comes from the condition for the self-consistency of the 
supersonic propeller model, see Eq. 7 in Ikhsanov 2002J. In 
this case the age of subsonic propellers turns out to be close 
to the typical ages of the SNRs with which they are presently 
associated. 

6. Discussion 

As shown above, some basic features of AXPs can be explained 
in terms of neutron stars in the state of subsonic propeller 
provided their mass and magnetic moment are M ~ 1 .4 Mq 
and fi ~ 2 x 10 30 Gcm 3 , and the strength of the wind is 
M c ~ 3 x 10 17 gs -1 . The required values of mass and dipole 
magnetic moment are typical for neutron stars. The required 
strength of the wind, however, is rather high and this condi- 
tion appears to be a strong restriction on the applicability of the 
proposed scenario. Indeed, for the wind strength ~ 10 17 gs~' 
to realize, the neutron star should be embedded in a medium of 
density 



10~ 16 Mn M\\ V 3 gcnT 3 . 



(13) 



Even in the case of a low relative velocity V re i ~ 10 6 cms _1 
(which is about the sound speed at the accretion radius of the 
star), the required number density of the surrounding gas, A^ e ~ 
10 s cirT 3 , remains significantly larger than the average density 
of material in SNRs and ISM. Hence, for the proposed scenario 
to be applicable one has to assume that AXPs are embedded in 
relatively dense regions. 

What are the physical conditions in these regions ? As fol- 
lows from X-ray observations, the column density in the ob- 
served AXPs is (0.4 - 2) x 10 22 cirT 2 (see e.g. Table 2 in 
Mereghetti et al. 2002). This indicates that the column density 
in the dense region is limited to 



N e L $ 4 x 10 21 cm" 2 , 



(14) 



where N e is the number density of the material in the dense 
region and L is the size of the region along the line of sight. 

If we require the life time of the accreting X-ray source to 
exceed the spin-down time scale of the subsonic propeller, one 
gets the following condition: 



, T - - ^acc^ns 

N e m p ® Z r d , 



(15) 



where <t> is the volume of the dense region. 

Finally, the total flux emitted by the material accumulated 
in the dense region should not exceed the fluxes detected from 
the infrared counterparts of AXPs. Under the conditions of in- 
terest, the dense regions are optically thin at the frequencies 
10 14 - 10 15 Hz. Hence, using the data presented in Fig. 2 of 
Israel et al. (2004i one finds: 



(16) 



W 2 r 1/2 <D ;S 6 x 10 58 , 



Combining Eq. (O - lO and setting T ~ 2000 K (the 
error of a factor of 2 in this parameter does not affect the final 
conclusions) we find that the subsonic propeller approach to the 
interpretation of AXPs is applicable if the corresponding X-ray 
sources are embedded in regions of volume <t> ~ 5 x 10 42 cm 3 , 
and average number density N e ~ 10 7 crrT 3 . The size of the 
regions along the line of sight is limited to L S 7 x 10 13 cm. 
Under these conditions the appearance of subsonic propellers 
is in agreement with properties of AXPs in both the X-rays and 
the infrared part of the spectrum. 

The question about the possible origin of such regions is 
rather complicated. In contrast to models built around the as- 
sumption that neutron stars in AXPs are in the state of accretor, 
the gravitational potential of the star within our scenario does 
not play a central role in the formation of its dense environ- 
ment. Indeed, the required size of the environment significantly 
exceeds the accretion radius of the star. This indicates that the 
formation of the dense region is rather independent of the prop- 
erties of the neutron star itself, but is related to the parameters 
of the supernova explosion or the properties of its progenitor or 
both. 

Three possible scenarios for the formation of dense envi- 
ronments around AXPs has recently been outlined by Marsden 



where T is the average temperature of the material in the dense 
region. 



et al. (2001): "fall-back" disk accretion, "pushed-back" disk 
accretion, and accretion involving high-velocity neutron stars. 
The fall-back disk accretion scenario within the subsonic pro- 
peller approach can be rejected since the material stored in the 
disk is assumed to be gravitationally bound to the star (this 
disk is situated within the accretion radius of the star, see e.g. 
Chatterjee et al. 2000 and references therein). 

The scenario in which a high-velocity neutron star captures 
material from the co-moving ejecta can be rejected on the same 
basis. Indeed, the number density in the expanding ejecta is 
n oc r 3 and therefore the life-times of AXPs within the sub- 
sonic propeller scenario turn out to be too short. The life-time 
could be longer if the material captured by the star were stored 
in a disk within its gravitational radius, but this assumption 
contradicts the basic statements of our model. 

Within the pushed-back disk scenario the dense environ- 
ment of the star is expected to form on a time scale of 10 3 yr. 
The total amount of material stored in such a disk is about 
0.4Mq (see Marsden et al. 120011 and references therein) 
and the disk is not gravitationally bound to the neutron star. 
Unfortunately, the basic properties of such a disk have been 
investigated rather poorly. For this reason, a detailed analysis 
of the applicability of this model to the subsonic propeller ap- 
proach is presently not possible. Nevertheless, the above men- 
tioned estimates of the mass and the time scale indicate that the 
pushed-back disk scenario could be a solution for the problem 
of the origin of the dense environment required in our model. 

It should also be noted that both the identification of the in- 
frared counterparts of AXPs (Israel et al. 2004) and the conclu- 
sion about the compact geometry of their host SNRs (Marsden 
et al. 2001} give some hints that there exists a dense environ- 
ment surrounding these objects. In the light of this, the require- 
ments of our model do not seem to be artificial but turn out to 
be rather reasonable from the observational point of view. 
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7. Conclusion 

The appearance of magnetized (B ~ 3 x 10 12 G) neutron stars 
in the state of subsonic propeller, which are situated in a strong 
(M c ~ 3 x 10 17 g s ) wind, is similar to that of AXPs in several 
important aspects. Namely, the rotational periods of these ob- 
jects are limited to the range of 5-15 s, and the spin-down rates 
are about 7 x 1CT 11 s s _1 . The emission of stars in this state is 
partly (5-10%) contributed to by the hot envelope surround- 
ing its magnetosphere and mainly (90-95%) by the process of 
accretion of material onto the stellar surface. The efficiency of 
the accretion process, which is governed by the reconnection 
of the magnetic field lines at the magnetospheric boundary, is 
high enough for the luminosity of the considered objects to be 
of order a few xlO 35 erg s . Finally, the estimated ages of neu- 
tron stars under the conditions of interest are close to 10 s yr. 

The above derived parameters allow us to suggest that the 
interpretation of AXPs within the subsonic propeller approach 
is rather promising for the future investigations. 
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